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The canonical front form Hamiltonian for non-Abelian SU(N) gauge theory 
in 3+1 dimensions is mapped non-perturbatively on an effective Hamiltonian 
which acts only in the Fock space of a quark and an antiquark. The ap¬ 
proach is based on the novel method of iterated resolvents and on discretized 
light-cone quantization, driven to the continuum limit. It is free of the usual 
Tamm-Dancoff truncations of the Fock space, rather the perturbative series are 
consistently resumed to all orders in the coupling constant. Emphasis is put on 
dealing with the many-body aspects of gauge field theory. Important is that 
the higher Fock-space amplitudes can be retrieved self-consistently from these 
solutions. 


1. THE STRUCTURE OF THE HAMILTONIAN 

In canonical field theory the four components of the energy-momentum vector 
P v commute and are constants of the motion. In the front form of Hamiltonian 
dynamics [jlj they are denoted p] by P u = (P + , P±, P~). Its spatial compo¬ 
nents P± and P + are independent of the interaction and diagonal in momen¬ 
tum representation. Their eigenvalues are the sums of the single particle mo¬ 
menta, P + = J^P + and Pl = Eal- Each single particle has four-momentum 
p M = (p + ,p_L,p~) and sits on its mass-shell p^p^ = m 2 . Each particle state 
u q” is then characterized by six quantum numbers q = (p + ,p±, A; c, /): the 
three spatial momenta, helicity, color and flavor. The temporal component 
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Fig. 1. - The Hamiltonian matrix for a SU(N)-meson and a harmonic resolution 
I\ = 4. Only vertex diagrams are included. Zero matrices are marked by a dot (■). 


P~ = 2 P + depends on the interaction and is a complicated non-diagonal oper¬ 
ator H[ |). It propagates the system in the light-cone time x + = x° + x 3 . The 
contraction of P M is the operator of invariant mass-squared, 

P"P M = P+P- - Pj 2 = H lc = H . (1) 

It is Lorentz invariant and refered to somewhat improperly but conveniently as 
the ‘light-cone Hamiltonian’ Plc> or shortly H. One seeks a representation in 
which H is diagonal P|'I'} = _E|\E'). Introductory texts ||, |[ |(| are available. 

The (light-cone) Hamiltonian is split for convenience into four terms: H = 
T+V+F+S. The kinetic energy T is diagonal in Fock-space representation, its 
eigenvalue is the free invariant mass squared of the particular Fock state. The 
vertex interaction V is the relativistic interaction per se. It has no diagonal 
matrix elements, is linear in g and changes the particle number by 1. The 
instantaneous interactions F and S are consequences of working in the light- 
cone gauge A + = 0. They are proportional to g 2 . By definition, the seagull 
interaction S conserves the particle number and the fork interaction F changes 
it by 2. As a consequence of working in light-cone representation and in the 
light cone gauge, the vacuum state has no matrix elemens with any of the Fock 
states: The vacuum does not fluctuate. 

The eigenvalue equation F[ jT) = P|\I/) stands for an infinite set of coupled 
integral equations ||. The eigenfunctions jT) are superpositions of all Fock- 
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space projections like (qq\\ k) or (qq g\^t). However, if one works with discretized 
light-cone quantization (DLCQ Q), one deals with a finite set of coupled matrix 
equations. As consequence of periodic boundary conditions the single particle 
momenta are discrete and the Fock states are denumerable and orthonormal. 
The Hamiltonian matrix is illustrated in Fig. |I] for the harmonic resolution 
K = P + L/(2n) = 4. Its rows and columns are the Fock states which are 
grouped into the Fock-space sectors qq, qqg,... , denumerated for simplicity 
by n = 1,2,..., 13. Within each sector one has many Fock states, which differ 
by the helicities, colors, and momenta of the single particles, subject to fixed 
total momenta P + and P±. Since P + has only positive eigenvalues and since 
each particle has a lowest possible value of p + , the number of particles in a Fock 
state is limited for any fixed value of K : The number of Fock space sectors is 
finite. The transversal momenta p± can take either sign. Their number must be 
regulated by some convenient cutoff, (Fock-space regularization, f.e. see ||). 
The most remarkable property of the Hamiltonian matrix is its block structure 
and its sparseness. Most of the blocks are zero matrices, marked by a dot (•) 
in the figure. Only those blocks in Fig. |l] denoted by the graphical symbol of 
a vertex are potentially non-zero, with the actual matrix elements tabulated 
elsewhere [j|. £j|. 

Here then is the problem, the bottle neck of any Hamiltonian approach in field 
theory: The dimension of the Hamiltonian matrix increases exponentially fast. 
Suppose, the regularization procedure allows for 10 discrete momentum states 
in each direction. A single particle has then about 10 3 degrees of freedom. A 
Fock-space sector with n particles has then roughly 10" _1 different Fock states. 
Sector 13 with its 8 particles has thus about 10 21 and sector 1 (qq) about 10 3 
Fock states. What one wishes instead of is to derive an effective interaction 
which acts only in the comparatively small qq- space, like a Coulomb interaction 
acts only between an electron and a positron. Loosely speaking, the aim of 
deriving an effective interaction can be understood as reducing the dimension 
in a matrix diagonalization problem from 10 21 to 10 3 ! 

The first attempt to formulate an effective interaction in field theory is the 
approach of Tamm and Dancoff (TDA Jt], ||). The present ‘method of iterated 
resolvents’ !, @ is closely related to TDA. Alternatively, one can apply a 
sequence of analytical and approximate unitary transformations in order to 
render the Hamiltonian matrix more and more ‘band diagonal’. Two methods 
have been proposed recently, the ‘similarity transform’ of Glazek and Wilson, 
see []ll|, and the ‘Hamiltonian flow equations’ of Wegner |l2|. All of these 
methods are under active research. Applications particularly by Trittmann et 
al [|~P~ij. |T1| . [l5|| , Brisudova et al |h], [TtJ, Jones et al |li| |n|, and Gubankova 
et al |2(| are presented elsewhere in these proceedings. Common to most of 
them is that one truncates the Fock space before approximation methods are 
applied. No truncation of the Fock space is needed with the method of iterated 
resolvents, as to be shown below. 
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2 . THE METHOD OF ITERATED RESOLVENTS 

The method of effective interactions is a well known tool in many-body physics 
which in field theory is known as the Tamm-Dancoff-approach ® !• The 
method is ideally suited for a field theory, because Fock-space sectors |z) appear 
in the most natural way. The Hamiltonian matrix can be understood as a 
matrix of block matrices, whose rows and columns are enumerated by i = 
1,2 .N like in Fig. |lj The eigenvalue equation can be written as a coupled 
set of block matrix equations: 

N 

^ {i\H\j) (?'|*F) = E (n|'F) for all* = 1, 2,..., N . (2) 

i=i 

In TDA, the rows and columns are split into a P-space, P = 11) (11, and the 
rest, the Q-space. Explicitly written out, the eigenvalue equation (||) becomes 
a set of two coupled block matrix equations: 

(P\H\P) (P\V) + (P\H\Q) {Q\V) = £<P|*>, (3) 

and (Q\H\P) (P\V) + (Q\H\Q) (Q\V) = E (Q |H>) . (4) 

Rewrite the second equation as (Q\E — H\Q) (QI’F) = ( Q\H\P) (P|H/), and 
observe that the quadratic matrix (Q\E — H\Q) could be inverted to express 
the Q-space wave-function (QI'F) in terms of (P^). But the eigenvalue E is 
unknown at this point. One introduces therefore a redundant parameter u>, 
and defines Gq(oj) = [{Q\u> — H IQ)]^ 1 . The so obtained effective interaction 

(P\H eS (u)\P) = (P\H\P) + (P\H\Q) G q (u) ( Q\H\P) . (5) 

acts only in the much smaller P-space: -ff e ff(w)|<i>fc(a;)) = Pfc(w) |*Ffc(w)). Vary¬ 
ing oj one generates a set of energy functions Ffc(w). All solutions of the fix- 
point equation Ek(co) = u are eigenvalues of the full Hamiltonian H. In fact 
one can find all eigenvalues of H despite the fact that the dimension of H e g is 
usually much smaller than that of the full H §0- The procedure is formal, 
since the inversion of the Q-space matrix is as complicated as its diagonaliza- 
tion. The adventage is that resolvents can be approximated systematically: 
The two resolvents 

W)= (Q\lu-T-U\Q) and G « )(w) = (Q\lo-T\Q) 

defined once with and once without the non-diagonal interaction U, are iden¬ 
tically related by Gq(u >) = Gq\u) + Gq\u>) U Gq(w), or by the infinite series 
of perturbation theory. Albeit exact in principle, the Tamm-Dancoff-approach 
(TDA) suffers from a practical aspect: The approach is useful only if one trun¬ 
cates the perturbative series to its very first term. This destroys Lorentz and 
gauge invariance, and requires a sufficiently small coupling constant. 
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Truncation can be avoided by ‘going backwards’. Reinterpreting the Q-space 
as the last sector TV (sector 13 in Fig. jlj) and the P -space with the rest, the 
above equations can be interpreted as to reduce the block matrix dimension 
from TV to TV — 1, with an effective interaction acting now in the smaller space 
of TV — 1 blocks. This procedure can be iterated, from TV — 1 to TV — 2, and so 
on, until one arrives at block matrix dimension 1 where the procedure stops: 
The effective interaction in the Fock-space sector with only one quark and one 
antiquark is defined unambiguously. One has to deal then with ‘resolvents of 
resolvents’, or with iterated resolvents [ jIo| . 

To be more explicit suppose one has arrived in the course of this reduction at 
block matrix dimension n, with 1 < n < TV. Denote the corresponding effective 
interaction H n {u). The eigenvalue problem reads then 

n 

^2{i\H n (uj)\j)(j\'i>{uj)) = E(u) , for i = 1,2,...,n. (7) 

.7=1 

Observe that i and j refer here to sector numbers, and that n refers to both, 
the last sector number, and the number of sectors. Now, like in the above, 
define the resolvent as the inverse of the Hamiltonian in the last sector 


{n\uj - H n {u)\n) 

n— 1 

G n MX>|ff n (w)U> (9) 

t=i 

The effective interaction in the (n — l)-space becomes then 

H n -i (w) = H n (u) + H n (cj)G n (u)H n (u;) (10) 

for every block matrix element (i\H n -i(ui)\j) . Everything proceeds like in 
above, including the fixed point equation E{ui) = uj. But one has achieved 
much more: Eq.(|lo|) is a recursion relation which holds for all 1 < n < TV! 
Since one has started from the bare Hamiltonian in the last sector, one has to 
convene that Hn = H. The rest is algebra and interpretation. 

Applying the method to the block matrix structure of QCD, as displayed 
in Fig. [I], is particularly easy and transparent. By definition, the last sector 
contains only the diagonal kinetic energy, thus = T 13 . Its resolvent is 
calculated trivially. Then H 12 can be constructed unambiguously, followed by 
Hu, and so on, until one arrives at sector 1. Grouping the so obtained results 
in a different order, one finds for the sectors with one gg-pair: 


H q ,= 

Hi 

= T-i+ vg 3 v + vg 3 vg 2 vg 3 v , 

( 11 ) 

Hqqg = 

h 3 

= T 3 + VG 6 V + VG 6 VG 5 VG 6 V + VG 4 V , 

( 12 ) 

Hqqgg = 

H 6 

= T 6 + VG W V + VG W VG 9 VG W V + VG 7 V. 

(13) 


G„(«) = 

thus (n | ^(w)) = 
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Fig. 2. — The three graphs of the effec¬ 
tive interaction in the qq- space. — The 
lower two graphs correspond to the chain 
U = VG 3 V, the upper corresponds to 
Ua = VG 3 VG 2 VG 3 V. Propagator boxes 
are represented by vertical dashed lines, 
with the subscript ‘n’ referring to the sec¬ 
tor numbers. 


The quark-gluon content of the sectors is added here for easier identification, 
in line with the notation in Fig. |I]. Correspondingly, one obtains for the sectors 
with two qq- pairs 

H qm = H 4 = T 4 + VG 7 V + VG 7 VG 6 VG 7 V , (14) 

H qqqqg = H 7 = Tr + VGnV + VGuVGwVGnV + VGsV. (15) 

In the pure gluon sectors, the structure is even simpler: 

Hgg = H 2 = T 2 + VG 3 V + VG 5 V , (16) 

Hggg= H 5 = T 5 + VG 6 V + VG 9 V . (17) 

Note that these relations are all exact. 

One is left with the eigenvalue problem in the gg-space, 

'52(q;q\H q g(uj)\q'-,q , )ip b (u)) = (lo) (q; q\ip b (a;)). (18) 

1' A' 

By construction, the eigenvalues are identical with those of the full Hamilto¬ 
nian and are enumerated by b = 1,2,.... They are the invariant mass 2 of a 
physical particle and its intrinsic excitations.The corresponding eigenfunctions 
(g;g|V’fc( w )) represent the normalized projections of |\I/) onto the Fock states 
|g; q) = b\d' q \vac). The effective Hamiltonian as given by Eq.( pd]) has two types 
of interactions which are illustrated diagrammatically in Fig. 2. In the first 
of them, in VG 3 V, the bare vertex interaction scatters the system virtually 
into the qq g-space, where it propagates under impact of the full interaction 
until a second vertex interaction scatters the system into the gg-space. The 
gluon can be absorbed on the same quark line as it was emitted, which con¬ 
tributes to the effective quark mass as illustrated by diagram U\^ 2 of Fig. ||. 
The gluon absorbed by the other quark as in diagram U\ 4 of Fig. providing a 
quark-antiquark interaction which cannot change the quark flavor. It includes 
all fine and hyperfine interactions, see for instance Jt 4|. The second term in 
Eq.dnJ), the annihilation interaction U a = VG 3 VG 2 VG 3 V, potentially pro¬ 
vides an interaction between different flavors and is illustrated in diagram t/yo 
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of Fig. |^. As a net result the interaction scatters a quark with helicity \ and 
four-momentum p into a state with \' q and p'. 

The knowledge of i/)b is sufficient to retrieve all desired Fock-space compo¬ 
nents of the total wave-function. The key is the upwards recursion relation 
Eq.(||). Obviously, one can express the higher Fock-space components (n| T) 
as functionals of %l> qq by a finite series of quadratures, i.e. of matrix multi¬ 
plications or of momentum-space integrations. One need not solve an other 
eigenvalue problem. This is quickly shown by way of example, by calculating 
the probability amplitude for a | gg)- or a |ggg)-state in an eigenstate of the 
full Hamiltonian. The first two equations of the recursive set in Eq.(^J) are 

<2|*> = G' 2 <2|7? 2 |1)<1|M>>, (19) 

and (3|’F) = G 3 (3|fT 3 |l)<l|*> + G 3 <3|i7 3 |2)<2|^). ( 20 ) 

The sector Hamiltonians H n have to be substituted from Eqs. dH) and (|l 6 |). In 
taking block matrix elements of them, the formal expressions are simplified con¬ 
siderably since most of the Hamiltonian blocks in Fig. |I| are zero matrices. One 
thus gets simply (2\H 2 \1) = {2\VG^V\1) and therefore (2|T) = G 2 VG^V(l\^). 
Substituting this into Eq. (|Ty) gives (3|>F) = G 3 V’(l|^ r ) + G 3 VG 2 VG^V{\\^). 
These findings can be summarized more succinctly as 

Hgg) = G gg VG qqg V \% q ), (21) 

and I ^qqg) = Gq q gV \ip qq ) + G qq g VG gg V G qq g V \4> qq ) ■ (22) 

The finite number of terms is in strong contrast to the infinite number of terms 

in perturbative series. Iterated resolvents sum the perturbative series to all 
orders in closed form. 

3. DISCUSSION AND PERSPECTIVES 

All of the above relations are exact and hold for an arbitrarily large K (Tc[| . 
They hold thus also in the continuum limit, where the resolvents are replaced 
by propagators and the eigenvalue problems become integral equations. The 
effective interaction is very simple in direct consequence of the structure of the 
gauge field Hamiltonian with its many zero matrices. No particle cut-off is 
required, and no assumption is made on the size of the coupling constant. The 
appoach is strictly non-perturbative. Instead of the inverting a huge matrix 
as in the Tamm-Dancoff approach, one has to invert only the comparatively 
small sector Hamiltonians H n . The instantaneous interactions can be included 
easily ex post by the rule that every intrinsic line in graph in Fig. |2j has to be 
suplemented with the corresponding instantaneous line p0| . Once the general 
structure of the effective interaction has been formulated, one can proceed 
with simplifying assumptions. For example, one can replace the full by the free 
propagators to get a selected set of perturbative diagrams. The major effort of 
the ongoing work | fll| is to find an approximation scheme which combines rigor 
with simplicity. Unfortunately the limited space prevents giving more details. 
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